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HIGHLIGHTS 


►  We  model  LiFeP04  battery  with  new  electrical  model. 

►  Parameters  in  battery  model  are  divided  into  loaded  and  unloaded  conditions. 

►  Transient  response  correction  is  proposed  to  formulate  the  parameters. 

►  We  validate  the  proposed  battery  model  with  experiment. 
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Lithium  ferro  phosphate  (LiFeP04)  batteries  are  stable  in  terms  of  their  thermal  and  chemical  charac¬ 
teristics  and  have  a  high  potential  for  wide  use  in  the  future.  A  battery  model  that  predicts  current 
—voltage  (/—  V)  performance  with  high  accuracy  is  vital  for  circuit  designers  so  they  can  design  circuits 
to  control  the  usage  of  the  battery  and  improve  its  runtime  and  its  safety  during  use.  In  this  paper, 
a  battery  model  with  two  resistance-capacitance  (RC)  parallel  networks  is  used  for  modelling  and  the 
values  for  parameters  are  divided  into  two  sets;  one  for  loaded  conditions  and  the  other  for  unloaded 
conditions.  A  methodology  of  transient  response  correction  is  proposed  for  formulating  the  parameters 
of  the  loaded  condition.  The  steps  for  parameter  extraction  are  discussed  and  the  proposed  battery 
model  is  validated  using  data  from  several  experiments.  The  comparison  between  experiment  and 
simulation  results  shows  that  the  proposed  model  is  capable  of  predicting  I—  V  performance  of  the 
battery  with  very  little  error. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Rechargeable  batteries  have  been  utilized  as  an  energy  storage 
element  for  portable  electronic  devices,  electric  vehicles  and 
distributed  generation  and  avionics  systems.  The  lithium-ion 
battery  has  been  the  dominant  type  used  in  portable  electronic 
devices  due  to  its  high  specific  energy  and  the  compact  size  [1  ].  The 
lithium  ferro  phosphate  (LiFeP04)  battery,  which  uses  phosphate  as 
its  cathode,  is  expected  to  be  widely  used  as  the  energy  storage 
element  in  the  future  due  to  its  advantages  in  thermal  and  chemical 
stability  [2]. 
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Even  though  battery  technology  has  significantly  improved  in 
recent  years,  modern  batteries  are  still  unable  to  meet  the  energy 
demands  of  current  portable  electronic  devices  [3,4].  Therefore, 
battery  management  is  important  for  battery-powered  devices 
because  it  reduces  power  dissipation  and  extends  the  operation 
time  of  the  battery.  In  order  to  have  an  effective  battery  manage¬ 
ment  system  for  a  LiFePCU  battery,  the  development  of  a  model  that 
can  accurately  predict  its  dynamic  behaviours  is  crucial  [5]. 

There  have  been  remarkable  efforts  committed  to  the  devel¬ 
opment  of  battery  models  over  the  past  decades.  The  three  cate¬ 
gories  of  battery  models  are  electrochemical  models,  mathematical 
models  and  electrical  models.  Because  the  electrochemical  models 
have  high  complexity  and  the  mathematical  models  are  unable  to 
provide  I—V  information,  these  models  are  not  well  suited  to  use  by 
circuit  designers  [3].  The  electrical  model  is  one  that  can  be  applied 
fairly  easily  by  circuit  designers  because  the  mathematical 
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Nomencalture 


List  of  symbols 

y  state-of-charge  of  battery,  % 

CN  nominal  capacity  of  battery,  Ah 

Rsr  series  resistance  for  full  charged  battery,  Q 
Ra  increment  of  series  resistance,  Q 
R i  and  R2  resistances  for  RC  parallel  networks  during  loaded 
conditions,  Q 

R\  and  R'2  resistances  for  RC  parallel  networks  during 
unloaded  conditions,  Q 
Vt  terminal  voltage  of  battery,  V 

Vki  and  Vk2  voltage  across  first  RC  parallel  networks  during 
loaded  conditions,  V 

Vkl  and  V’k2  voltage  across  first  RC  parallel  networks  during 
unloaded  conditions,  V 

iki  and  ik2  time  constants  of  the  transient  response  during 
loaded  conditions,  s 

T'kl  and  t[<2  time  constants  of  the  transient  response  during 
unloaded  conditions,  s 

A  and  C  parameters  for  unloaded  conditions,  V 
B  and  D  parameters  for  unloaded  conditions,  s_1 
E  and  G  parameters  for  unloaded  conditions,  V 
Fand  H  parameters  for  unloaded  conditions,  s_1 
M  and  N  gains  for  the  gain  correction 
L  gain  for  the  time  constant  correction 

Tsi  and  Ts2  settling  times  of  the  transient  response,  s 


equations  can  be  implemented  and  the  parameters  derived  using 
electrical  circuits  [6].  Furthermore,  the  I—V  information  can  be 
readily  extracted  from  the  electrical  model. 

The  most  simple  electrical  battery  model  is  a  circuit  with 
a  controlled  voltage  source  in  series  with  a  resistor,  however,  due  to 
the  transient  response  characteristic  associated  with  the  electro¬ 
chemical  process  of  a  rechargeable  battery,  a  second  type  of  elec¬ 
trical  battery  model  called  the  voltage-resistance-capacitance 
(VRC)  model  has  been  developed.  The  VRC  model  is  a  circuit  with 
a  controlled  voltage  source  in  series  with  a  resistor  and  a  RC  parallel 
network  as  used  in  [6-12].  The  VRC  model  is  improved  by  utilizing 
two  RC  parallel  networks  as  proposed  by  Chen  and  Mora  [3].  It  is 
useful  in  characterizing  both  the  short  and  long  time  transient 
responses  of  the  electrochemical  process.  Circuit  diagram  of  the 
VRC  model  with  two  RC  parallel  networks  is  shown  in  Fig.  1.  The 
VRC  model  with  two  RC  parallel  networks  has  been  shown  to  have 
achieved  high  accuracy  [3,4,13,14]. 

The  parameters  for  the  VRC  model  have  been  identified  based  on 
the  voltage  curve  for  the  unloaded  condition.  There  has  been 
a  disagreement  about  parameter  identification  and  Norian  has 
pointed  out  that  the  parametric  values  for  loaded,  unloaded  and 
charged  conditions  are  not  necessarily  the  same  [8].  The  parameters 
obtained  from  the  unloaded  conditions  are  not  suitable  for  the 
voltage  response  under  loaded  conditions.  Flence,  the  transient 
response  correction  should  be  added  to  the  existing  electrical  battery 
model  so  that  it  can  be  used  for  LiFePC^  batteries.  The  transient 
response  correction  technique  is  straightforward  and  therefore 
simplifies  the  process  of  parameter  extraction.  The  technique  can  be 
used  to  improve  the  accuracy  of  the  existing  battery  model. 

In  this  paper,  a  battery  model  with  different  parametric  values  of 
RC  parallel  networks,  for  loaded  and  unloaded  conditions,  is 
developed  using  transient  response  correction.  The  quality  of  the 
developed  battery  model  is  validated  using  data  obtained  from 
actual  battery. 


Fig.  1.  Electrical  battery  model  which  proposed  by  Chen  and  Mora  [3]. 


2.  Improved  model 

Fig.  2  illustrates  the  proposed  electrical  battery  model.  It  uses 
the  same  electrical  circuit  as  the  VRC  model  but  with  different 
parametric  values  of  RC  parallel  networks  for  loaded  and  unloaded 
conditions.  The  terminal  voltage  and  loaded  current  are  repre¬ 
sented  by  Vt  and  /L  respectively.  The  open-circuit  voltage  (OCV)  is 
represented  by  a  voltage  source  which  is  dependant  on  the  state-of- 
charge  (SoC)  of  the  battery.  SoC  represents  the  remaining  capacity 
of  the  battery  as  denoted  by  Eq.  (1 )  where  y{t\f)  is  the  SoC  at  time  tk, 
y(to)  is  the  initial  SoC  and  Cn  is  nominal  capacity  of  the  battery. 

tk 

Y(tfc)  =  r(to)  -  /  k(t)dt  ,  to  <  t  <  tk  (1) 

to 

The  series  resistance  of  the  battery  is  characterized  by  both  Rsr 
and  Ra.  Rsr  is  used  to  represent  the  series  resistance  when  the 
battery  is  fully  charged  and  its  value  is  constant,  however,  the  series 
resistance  is  found  to  be  larger  when  the  battery  is  loaded  under 
partially  charged  conditions.  Flence,  a  resistor  Ra  is  used  to  repre¬ 
sent  the  increment  in  the  series  resistance.  In  this  formulation,  the 
value  of  Ra  is  dependant  on  the  initial  SoC  when  the  battery  is 
loaded  and  the  value  of  Ra  is  zero  when  the  battery  is  fully  charged 
(SoC=100%). 

The  transient  characteristic  of  the  battery’s  electrochemical 
process  is  characterized  by  using  two  RC  parallel  networks.  For 
loaded  conditions,  the  parameters  of  RC  parallel  networks  are  R\, 
R2 ,  Ci  and  C2  as  shown  in  Fig.  2(a)  whereas  the  parameters  of  RC 
parallel  networks  are,  R\,  R'2 ,  C\  and  C2  for  unloaded  conditions  as 
shown  in  Fig.  2(b).  The  details  on  parameter  identification  are 
discussed  in  Section  4. 

3.  Experiment  set  up  and  procedures 

3.1.  Experiment  set  up 

The  block  diagram  and  the  experiment  set-up  for  the  battery 
test  are  shown  in  Figs.  3  and  4.  The  test  system  consists  of  a  LiFePCU 


Fig.  2.  Proposed  electrical  battery  model  for  (a)  loaded  conditions  and  (b)  unloaded 
conditions. 
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Fig.  3.  Block  diagram  of  battery  tests. 


Fig.  4.  Equipment  set-up  of  battery  tests. 

battery,  whose  nominal  capacity  and  nominal  voltage  are  18  Ah  and 
3.2  V  respectively.  The  DAQ  NI9219  which  produced  by  National 
Instruments  is  interfaced  with  LabVIEW  2010  and  a  computer  to 
collect  data  from  the  experiment.  The  battery  tests,  which  require 
a  constant  current  load,  are  carried  out  by  using  electronic  load 
IT8514C  (produced  by  ITECH)  with  ratings  of  120  V,  240  A  and 
1200  W. 

3.2.  Battery  test  procedure 

Pulse  discharge  (PDT)  and  continuous  discharge  (CDT)  tests  are 
performed  on  the  LiFePCU  battery  for  parameter  identification  and 
model  validation.  The  PDT  test  consists  of  a  sequence  of  constant 
current  discharge  and  rest  periods.  The  battery  is  loaded  with  9  A 
(0.5  C)  of  current  for  720  s  to  reduce  10%  of  the  battery’s  SoC.  After  4  h 
of  rest,  the  same  discharge  pulse  is  applied  until  the  battery  voltage 
reaches  2.0  V.  In  this  configuration,  the  battery  is  unloaded  for  a  long 


|  ^—Battery  Terminal  Voltage  [ 


0  0.5  1  1.5  2  2.5 


Fig.  5.  Voltage  response  and  current  profile  for  PDT  test. 


Time  (s) 

Fig.  6.  Voltage  response  and  current  profile  for  CDT  test. 

period  (10  times  longer  than  the  RC  transient)  to  accurately  deter¬ 
mine  the  parameters  of  the  battery  model  as  suggested  in  [15].  In 
this  modelling  technique,  battery  is  required  to  be  in  full  relaxed 
condition,  thus  highly  transient  current  is  not  suitable  for  parameter 
identification.  Also,  in  order  to  avoid  permanent  damage  caused  by 
over  usage,  the  battery  voltage  is  strictly  controlled  to  be  higher  than 
2.0  V.  Several  PDT  tests  are  carried  out  with  different  current  pulses 
in  order  to  investigate  the  parameters  for  several  current  loads.  On 
the  other  hand,  the  CDT  test  is  done  by  discharging  the  battery  with 
continuous  constant  current.  Continuous  discharge  tests  for  9  A 
(0.5  C),  13.5  A  (0.75  C)  and  18  A  (1  C)  of  current  are  made.  The  typical 
voltage  response  and  the  current  profile  for  PDT  and  CDT  battery 
tests  are  shown  in  Figs.  5  and  6  respectively. 

4.  Model  extraction  and  transient  response  correction 

In  the  developed  model,  all  the  parameters  are  expressed  as 
a  function  of  SoC.  Because  the  value  of  SoC  is  decreasing  during 
loaded  conditions,  the  identification  of  parametric  equations  is 
difficult.  Therefore,  the  model  extraction  is  initiated  from  the 
unloaded  conditions  with  a  constant  value  of  SoC.  After  the  values 
of  parameters  are  determined  from  the  unloaded  conditions, 
a  transient  response  correction  is  employed  in  order  to  determine 
the  value  of  parameters  for  the  loaded  conditions. 

4.1.  Parameter  identification  for  the  unloaded  conditions 

As  seen  in  Fig.  2(a),  the  value  of  loaded  current,  4,  is  zero  for  the 
unloaded  conditions.  Consequently,  the  values  of  series  resistances 
Rsr  and  Ra  are  not  important  for  the  unloaded  conditions.  The 
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Fig.  8.  Relaxation  curve  at  80%  SoC  for  9  A  (0.5  C),  13.5  A  (0.75  C)  and  18  A  (1  C)  PDT 
tests. 


4.1.1.  Open-circuit  voltage,  OCV 

OCV  for  a  certain  SoC  is  determined  by  the  terminal  voltage  of 
the  battery  after  the  battery  is  unloaded  for  a  long  period.  As  sug¬ 
gested  in  [15],  the  OCV  can  be  precisely  determined  after  the 
battery  is  unloaded  for  at  least  10  times  the  time  of  the  RC  transient 
response.  In  this  paper,  the  OCV  is  determined  after  the  battery  is  at 
rest  for  4  h.  The  OCV— SoC  curve  is  shown  in  Fig.  7  and  the  rela¬ 
tionship  is  expressed  by  Eq.  (3),  where  y  represents  SoC. 

OCV(y)  =  (-3.001  x  10-n)r6+  (9.85  x  lCT9)?5 

-  (l.291  x  lCT6)?4  +  (8.763  x  lCT5)?3 

-  0.003325?2  +  0.07002?  +  2.607  (3) 


terminal  voltage,  Vt,  is  denoted  as  Eq.  (2).  Vkl  and  V[a  are  the 
voltages  across  the  RC  parallel  networks  during  the  unloaded 
condition,  whereas  VkiM  and  V\q[ol)  are  the  voltage  levels  across 
the  RC  parallel  networks  at  the  beginning  of  unloaded  condition,  a. 
T'kl  and  t[<2  represent  the  time  constants  of  the  transient  response  of 
Vkl  and  V[a  respectively.  In  this  paper,  the  parameter  identification 
is  done  using  MATLAB  curve-fitting  tools. 


v't(t)  =  ocvhr(t)]-  ^(0-^(1) 
V*lW -M«)exp(=^)  t>a 

K 2(0  =  via(a)exp^{t^a)\ 


(2) 


4.1.2.  Parameters  of  transient  response  for  the  unloaded  conditions 
The  voltage  response  for  the  unloaded  condition  is  a  combina¬ 
tion  of  the  transient  responses  having  both  a  short  time  and  long 
time  constant.  In  this  paper,  two  RC  parallel  networks  are  used  to 
represent  the  voltage  responses.  In  order  to  simplify  the  modelling 
method,  the  equations  for  Vkl  and  V[a  can  be  expressed  as  Eq.  (4), 
where  A,  B,  C,  and  D  represent  the  desired  values  of  Vki(a),  i]^1, 
Vk2 (a)  and  respectively.  Thus,  the  parameter  identification  of 

R\,  R'2,  q  and  C2  can  be  avoided.  It  is  important  to  note  that  the 
parameters  A,  B,  C  and  D  are  independent  of  discharge  current,  4, 
since  the  relaxation  curves,  as  shown  in  Fig.  8,  are  similar  for  all 
current  rates. 


KAO  =  Aexp(—B(t  -  a)) 
V^2(t)  =  Cexp(-D(t  -  a))  ’ l> 


(4) 


where, 


where, 

^  =  R\  X  =  R'2xC'2 


A  =  Vkl(a),C  =  Vk2(a) 

1  ~  1 


B  = 


,D  = 


kl 


1<2 


> 

o 


Fig.  9.  Extracted  parameters  for  unloaded  conditions. 
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Fig.  10.  Ra- SoC  relationship. 


This  trend  indicates  that  the  process  of  parameter  extraction  for 
A ,  B ,  C  and  D  is  straightforward.  The  curve  of  parameters  versus  SoC 
is  shown  in  Fig.  9.  The  identified  parameters  can  be  expressed  by 
Eq.  (5).  The  goodness-of-fit  statistics  of  the  curve  fitting  is  tabulated 
in  Table  1. 


4.2.1.  Series  resistances,  Rsr  and  Ra 

The  series  resistance  is  represented  by  both  Rsr  and  Ra  where  Rsr 
represents  the  series  resistance  for  a  fully-charged  battery  and  Ra 
represents  the  increment  of  series  resistance  for  other  SoC.  Series 
resistance  can  be  calculated  from  the  instantaneous  voltage  drop 
when  the  load  is  applied.  For  this  paper,  the  value  of  Rsr  is 
0.005889  Q  while  the  Ra  is  varied  with  SoC  as  shown  in  Fig.  10.  The 
average  value  of  Ra  can  be  expressed  by  Eq.  (7). 


J?a[y]  =  0.1457e~°  079367  (7) 

4.2.2.  Parameters  of  transient  response  for  the  loaded  conditions 
The  terminal  voltage,  Vh  of  the  battery  for  the  loaded  conditions 
can  be  expressed  by  equation  (6).  During  the  loaded  conditions,  the 
SoC  of  the  battery,  7,  will  decrease.  Because  all  the  parameters  of 
transient  response  are  dependant  on  the  SoC,  the  functions  for 
Vki(y).  142(7).  Tki (t)  and  Tk2(r)  cannot  be  accurately  determined 
from  the  voltage  curve.  The  transient  response  correction  must  be 
applied  to  determine  the  value  of  the  parameters.  An  alternative 


A(y)  =  -(2.5348  x  10-8)y4  +  (5.8161  x  KT6)?3  -  (4.6837  x  lO”4)?2  +  (1.5226  x  \Q-2)y  -  0.13906 
B(y)  =  -(5.6107  x  10-6)y  +  (8.1434  x  lCT4) 

C(y)  =  (3.9383  x  10-8)y4  -  (8.3679  x  10-6)y3  +  (6.4756  x  10~4)y2  -  (2.212  x  10-2)y  +  0.32153 
D(y)  =  -(3.6536  x  10-5)y  + 0.012037 


4.2.  Parameter  identification  for  the  loaded  conditions 

The  voltage  response  of  the  loaded  conditions  consists  of  an 
instantaneous  voltage  drop  and  a  series  of  transient  responses.  In 
this  context,  because  the  loaded  current  is  not  zero,  the  series 
resistances,  Rsr  and  Ra  must  be  determined.  The  terminal  voltage,  Vt, 
can  be  expressed  by  Eq.  (6).  V\ ci  and  V\a  are  the  voltages  across  the 
RC  parallel  networks  during  the  loaded  condition  while  ti<i  and  ti<2 
represent  the  time  constants  of  the  transient  response  of  V\ a  and 
142  respectively.  The  parameters  of  the  transient  responses  are  then 
identified  using  a  transient  response  correction  technique. 


way  to  determine  parameters  for  loaded  conditions  is  optimisation 
technique  as  discussed  in  [16-18].  Generic  algorithm  is  applied  in 
optimisation  technique  to  generate  parameter  set  for  loaded 
condition,  however,  in  transient  response  correction,  the  parameter 
identification  is  straightforward.  The  equations  for  V\ a  and  V\a  are 
simplified  as  Eq.  (8). 


vM  (t) 
Vk2(t) 


E(  1  -  exp(  -  F(t  -  t0))) 
G(1  -  exp(  -  H(t  -  t0))) 


,t0  <t  <  a 


where, 


(8) 


Vt (t)  =  OCV[Y(t)]  -  Vkl(t)  -  Vk2(t)-Ih(t)  X  (Rsr  +Ra) 
where, 

Vki(t)  =  Ii(t)  xR,(l  -exp(^M)) 

V  \  kl  /  /  +r\  <C  t  <C  (Y 

vE(()  =  «Ox«2(i-exr(^))'0 

and, 

T1<1  =  ^1  X  Q5t1<2  =  ^2  x  C2 


(6) 


Table  1 

Goodness-of-fit  of  curve  fitting. 


Goodness-of-fit 

A 

B 

C 

D 

SSE 

3.758xl0-6 

1.7xl0-7 

2.688xl0-6 

2.551  xlO-5 

R-square 

0.9811 

0.3415 

0.9984 

0.1278 

Adjusted 

0.9432 

0.2098 

0.9951 

-0.04,665 

R-square 

RMSE 

0.001371 

0.0001844 

0.001159 

0.002259 

Fig.  11.  Flow  chart  of  transient  response  correction. 
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Fig.  12.  Mismatch  between  experiment  data  and  simulation  result  when  the  setting  of 
E  =  A,  F  =  B.  G  =  C,  and  H  =  D. 


E  =  k(t)  x  R,,G  =  k(t)  x  R2 

F  -  1  ,H  .  1 

Tkl  t1<2 

4.2.3.  Transient  response  correction  method 

The  flow  chart  of  the  transient  response  correction  method  is 
shown  in  Fig.  11.  The  steps  of  the  transient  response  correction 
method  are  as  described  as  below. 

4.2.3.  t.  Gain  correction.  The  gain  correction  is  started  by  assuming 
that  the  transient  response  parameters  for  the  loaded  conditions  E 
to  H  have  the  same  value  as  A  to  D.  Because  the  values  of  Vki(to)  and 
V\a{to)  for  the  loaded  condition  are  unknown,  they  are  set  to  zero 
whilst  the  values  of  14i(a)  and  V\a{a)  are  set  to  A  and  C, 


to  C( 

Time(s) 

Fig.  13.  Voltage  curve  of  Vki  and  Vk2  (a)  before  and  (b)  after  gain  correction. 


Fig.  14.  Comparison  between  experiment  data  and  simulation  result  after  gain 
correction,  with  the  setting  of  E  =  M  x  A,  F  =  B,  G  =  N  x  C  and  H  =  D. 


respectively,  for  the  unloaded  conditions.  An  example  of  mismatch 
between  experiment  and  simulation  results  is  shown  in  Fig.  12.  The 
mismatch  between  experiment  and  simulation  results  proves  that 
the  parameters  of  transient  response  for  loaded  and  unloaded 
conditions  are  not  the  same.  Therefore,  a  transient  response 
correction  method  is  needed. 

In  order  to  perform  a  transient  response  correction,  the  voltage 
across  the  RC  parallel  networks  is  measured.  Voltage  responses  for 
the  voltage  across  the  RC  parallel  networks  before  and  after  the 
transient  response  correction  are  shown  in  Fig.  13.  Because  the 
voltage  curves  of  and  V^2  represent  the  desired  voltage  curve 
for  the  voltage  across  the  RC  parallel  networks  in  the  relaxation 
condition,  the  Vki  and  V\a  curves  should  match  the  curves  of  Vj^ 
and  V[a  at  time  a ,  however,  as  can  be  seen  in  Fig.  13(a),  a  voltage 
gap  between  I4i  and  V[a  is  found.  In  order  to  meet  the  desired 
value  of  “A”  during  the  unloaded  condition,  the  value  of  Vki(cO 
should  be  increased  to  the  value  of  “A”.  The  correction  can  be  made 
by  multiplying  the  value  of  E  with  a  gain,  M,  as  shown  in  Eq.  (9).  In 
Fig.  13(a),  there  is  no  voltage  gap  for  V\a,  however,  the  same 
concept  is  applied  by  multiplying  G  with  the  gain  of  N  if  a  voltage 
gap  exists  between  V\a  and  V[a.  The  value  of  Ri  and  R2  are  deter¬ 
mined  by  using  equation  (10)  after  the  value  of  E  and  G  are  updated 
to  Enew  and  Gnew  respectively.  The  mismatch  between  experi¬ 
mental  and  simulation  results  can  be  effectively  solved  as  shown  in 
Fig.  14. 


Fig.  15.  Example  of  mismatch  between  experiment  data  and  simulation  result  after  the 
gain  correction. 
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Fig.  16.  Comparison  between  experiment  data  and  simulation  result  (a)  before  and  (b) 
after  time  constant  correction,  with  the  setting  of  E  =  A,F  =  LxD,G  =  C  and  H  =  L  x  D. 


Enew  —  1W  X  £  —  M  x  IK 

Gnew  =  NxG  =  NxC 
Therefore, 


Vki(«) 

Vki(«) 

Via  («) 


xA 

x  C 


k{t)  xRi 

k(C)  X  R2 


ri 


Ri 


Enew 

w 

Gnew 

w 


(9) 


(10) 


42.3.2.  Time  constant  correction.  Example  of  the  mismatch 
between  experiment  and  simulation  results  after  the  gain  correc¬ 
tions  is  shown  in  Fig.  15.  In  this  paper,  this  type  of  mismatch  can  be 
found  at  the  high  SoC  region  (SoC  more  than  60%).  As  can  be  seen, 
although  the  voltage  curve  of  the  experiment  data  and  simulation 
result  meet  at  the  same  point  at  time  a,  the  rate  of  voltage  drop  is 
different.  The  rate  of  voltage  drop  for  the  experiment  data  is  faster 
than  the  simulation  result  implying  a  deviation  in  the  time 
constant.  Hence,  a  time  constant  correction  is  made  by  shortening 
the  time  constant  value  of  the  transient  response. 


Table  2 

The  average  value  of  correction  gains  in  several  SoC. 


SoC  (%) 

M 

N 

L 

90  <  y  <  100 

3.1667 

80  <  7  <  90 

1.6 

70  <  7  <  80 

2.4 

60  <  7  <  70 

1.5 

50  <  7  <  60 

3.51,693 

1.01,228 

40  <  7  <  50 

2.9898 

1.036 

7  <  40 

2.1 

1.00,944 

Fig.  17.  Comparison  between  simulation  result  with  the  experiment  data  of  PDT  test 
for  (a)  9  A  (0.5  C),  (b)  13.5  A  (0.75  C)  and  (c)  18  A  (1  C). 

The  time  constant  correction  is  started  by  removing  the  gain 
corrections,  M  and  N.  Then,  the  time  constant  correction  is  made  by 
setting  both  iki  and  z\a  to  be  equal  to  the  value  of  the  shorter  time 
constant  in  the  unloaded  condition,  t{<2.  The  values  of  Fand  H  are 
set  to  be  equal  to  D  and  the  transient  response  of  the  battery  has 
now  become  first  order. 

Processes  of  time  constant  correction  are  illustrated  in  Fig.  16. 
After  the  time  constant  is  set  to  D-1,  the  comparison  between 
experiment  and  simulation  results  is  made  as  shown  in  Fig.  16(a). 
The  comparison  shows  that  the  settling  time  of  the  transient 
response  is  different.  The  correction  of  the  time  constant  is  made  by 
using  a  time  constant  gain,  L,  as  expressed  in  Eq.  (10).  The  result  of 
the  comparison  after  the  time  constant  correction  is  shown  in 
Fig.  16(b).  The  result  shows  that  the  correction  effectively  solves  the 
problem  of  time  constant  mismatch. 


Table  3 

RMS  modelling  errors  in  PDT  tests. 

Current  (A) 

RMS  error  (V) 

9 

0.0035 

13.5 

0.0116 

18 

0.0079 
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Fig.  18.  Comparison  between  simulation  result  and  experiment  data  of  CDT  test  for  (a) 
9  A  (0.5  C),  (b)  13.5  A  (0.75  C)  and  (c)  18  A  (1  C). 

Average  values  of  gain  correction  in  several  SoC  for  the  loaded 
current  between  9  A  (0.5  C)  to  18  A  (1  C)  is  tabulated  in  Table  2.  For 
this  case,  it  was  found  that  the  time  constant  correction  is  imper¬ 
ative  for  the  SoC  region  higher  than  60%  whereas  the  gain  correc¬ 
tion  is  applied  for  the  SoC  region  of  60%  and  below. 

5.  Model  validation 

In  order  to  validate  the  developed  model,  the  experiment  data  of 
PDT  and  CDT  tests  are  compared  with  the  simulation  results. 
Comparison  between  experiment  data  and  with  the  simulation 
result  of  the  PDT  tests  is  shown  in  Fig.  17.  The  load  profiles  are: 


Table  4 

RMS  modelling  errors  in  CDT  tests. 


Current  (A) 

RMS  error  (V) 

9 

0.0094 

13.5 

0.0124 

18 

0.0125 

Fig.  19.  Current  profile  in  random  pulse  test. 


(a)  9  A  (0.5  C)  constant  current  pulse  with  a  720  s  duty  cycle  and 
a  period  of  1  h  per  cycle.  (10%  of  SoC  reduction  per  cycle). 

(b)  13.5  A  (0.75  C)  constant  current  pulse  with  a  480  s  duty  cycle 
and  a  period  of  30  min  per  cycle.  (10%  of  SoC  reduction  per 
cycle). 

(c)  18  A  (1  C)  constant  current  pulse  with  a  720  s  duty  cycle  and 
a  period  of  1  h  per  cycle.  (20%  of  SoC  reduction  per  cycle). 

There  is  a  close  agreement  between  experiment  and  simulation 
results  as  shown  in  Fig.  17.  The  0.03  V  of  maximum  error  for  the  PDT 
tests  is  found  at  the  beginning  of  discharging  (within  a  short  time). 
As  discussed  in  Section  4.2.3.2,  the  RC  transient  responses  are 
adjusted  to  be  first  order  when  the  time  constant  correction  is 
made.  Therefore,  the  simulation  result  is  slightly  different  from  the 
experiment  data  due  to  the  small  deviation  between  first  order  and 
second  order  transient  responses.  The  RMS  modelling  errors  in  PDT 
tests  are  tabulated  in  Table  3. 

On  the  other  hand,  the  comparison  between  experiment  and 
simulation  results  for  several  CDT  test  are  shown  in  Fig.  18.  The 
comparative  analysis  shows  that  experiment  and  simulation  results 
are  well  matched.  The  maximum  error,  which  happens  at  the 
beginning  of  discharging,  is  0.018  V  and  is  caused  by  the  small 
deviation  between  first  order  and  second  order  transient  responses. 
The  RMS  modelling  errors  in  CDT  tests  are  tabulated  in  Table  4. 

The  developed  model  is  further  validated  with  random  pulse 
test.  In  this  test,  battery  is  loaded  with  various  current  rates  as 
shown  in  Fig.  19,  which  including  3  A  (0.1667  C),  6  A  (0.3333  C),  12  A 
(0.6667  C),  and  36  A  (2  C)  of  current.  Moreover,  some  of  the 
charging  conditions  are  also  included  in  the  test.  The  comparison 
between  experiment  and  simulation  results  for  random  pulse  test 
is  shown  in  Fig.  20.  The  comparative  analysis  shows  that  experi¬ 
ment  and  simulation  results  are  well  matched  with  RMS  modelling 
error  of  0.0189  V. 


Fig.  20.  Comparison  between  simulation  result  and  experiment  data  of  random  pulse 
test. 
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Based  on  the  good  match  between  experiment  and  simulation 
results,  the  developed  model  is  validated.  This  proves  that  the 
transient  response  correction  technique  is  capable  of  reducing  the 
mismatch  between  the  model  and  the  experiment  data. 

6.  Conclusions 

In  this  paper,  an  accurate  and  simple  electrical  battery  model  is 
successfully  developed.  The  modification,  to  include  the  transient 
response,  of  an  existing  electrical  battery  model  is  successfully 
applied  to  reduce  the  modelling  error  in  both  loaded  and  unloaded 
conditions.  The  close  agreement  between  simulation  and  experi¬ 
ment  results  shows  that  the  characteristics  of  a  LiFePCU  battery  can 
be  accurately  simulated  with  the  developed  model.  Moreover,  the 
simple  approaches  in  model  extraction  enable  the  researcher  to 
easily  build  up  a  higher  accuracy  battery  model.  With  the  proposed 
model,  the  battery  modelling  method  can  be  simplified  and  will 
eventually  help  in  the  development  of  battery  management  systems 
and  energy  storage  systems  for  green  technology  applications. 
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